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The intensity of the overlap of a quantum state with all its phase space translations defines 
its quantum correlations. In the case of pure states, these are invariant with respect to Fourier 
transformation. The overlaps themselves are here studied in terms of the Wigner function and its 
Fourier transform, i.e., the characteristic function or chord function. Unlike the Wigner function, 
the chord function need not be real, but eventual symmetry with respect to reflections about a phase 
' space point may relate these representations. Semiclassical approximations for the "classical-like" 

region of small chords and for large chords are derived. These lead to an interpretation of the Fourier 
invariance in terms of conjugate chords. The interrelation of large and small (sub-Planck) phase 
space structures previously noted in the literature are thus reinterpreted. 
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I. INTRODUCTION 

> 

. Phase space correlations in quantum states have quite different properties than correlations in classical distribu- 

' tions functions. At the quantum level, the uncertainty principle imposes, on one hand, limitations to the possible 

i: distribution functions and at the same time creates strong relations between small and large scales that result in 

I ■ peculiar properties. Thus, for example, it has been realized that the phase space area A over which a pure state 

, extends determines the minimal size of the high frequency oscillation structures SA in the Wigner function by a kind 

' of complementarity relationship. 



A-SA>{2TTnf (1) 



(for one degree of freedom), thus characterizing these structures as "sub-Planck". This relationship between possibly 
^ macroscopic areas A and sub-Planck areas 6 A is entirely due to the finiteness of Planck's quantum of action and is 
^ not present when considering the correlations of classical distribution functions. In general, and especially for systems 
with more than one degree of freedom, we shall show that it is more revealing to relate a given large displacement to 
^ ' a specific small scale oscillation, rather than to relate areas. 
• I— I . This complementarity between small and large scales in quantum distribution functions can be accessed by the 
' Wigner function and its Fourier transform. Moreover, the latter, sometimes referred to as the characteristic function, 
^ , the generating function, or simply the chord function, also lies at the core of a full represention of quantum mechanics, 
- - ' on a par with the Weyl- Wigner representation itself. Further considerations about the corresponding conjugate 
classical phase spaces as well as the notation is found in Appendix ^ In this paper we compare the way that 
both these representations exhibit quantum coherences on all scales, that are overlaid on purely classical structures. 
The copious previous litterature on the Wigner functions allows us to concentrate on the properties of the chord 
representation. 

In Section^ we present the intrinsic definition of phase space correlations for a general density operator p to be 
studied. These are then related to the chord function and the Wigner function in Section Hill Then in Section HVI 
we exploit the general interrelation of these functions in the case of pure states and present some simple examples. 
Section is dedicated to the theory of small chords. This is the classical region, corresponding to a caustic of the 
Wigner function. The simple approximation thus obtained is rederived in the next section within the full semiclassical 
theory valid for large chords. Squaring the semiclassical chord function in section rVIII we interpret the invariance of 
the correlations of a pure state with respect to the Fourier transform: This results from a geometrical conjugacy of 
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chords defined on the classical phase space structures. Thus the complementarity of large and small chords arises as 
a limit among finite conjugate chords. The discussion in the last section recapitulates the full picture for extended 
pure states including "ergodic states" of chaotic Hamiltonians. 



II. PHASE SPACE CORRELATIONS 

The correlation between two states of a quantum system, represented by their density operators, pA and ps, can 
be given an intrinsic definition as 

C^^ ^ JlM^ . (2) 
^/trp\ tr/5| 

By the Schwartz inequality this quantity is always less than unity and for pure states, p — \'ijj){'ip\, it reduces to 
Cab = |('0a|V'b)P- When the pair of states are obtained by the unitary evolution of slightly different hamiltonians 
from the same initial state, the behaviour in time of Cab reflects the regular or chaotic nature of the underlying 
classical motion [^|^^. If the states are related by a unitary transformation generated by an Hermitian operator K, 

with a a real parameter, the small-a behaviour of the correlation is easy to ascertain: 

^. ,trpe-'"^^''pe"^^^^ tT\p^ - pK pK] , a^tr\p,K]'^ ^ 

C[a) = — w 1 7T = 1 H 7 — < 1 . (4) 

If p represents a pure state, this quadratic behaviour relates C{a) to the dispersion of the generator K, i.e., {K^) — {K)"^ 

In what follows we will be concerned with the phase space correlations of a quantum state produced by unitary 
translations in the L-dimensional phase space (p, q) = {pi, . . . gi, . . . , g^). Denoting the corresponding quantum 
operators by (p, q) , the translation operators read 



— exp 



(5) 



where the vector, or chord, ^ = (CpiCg) represents an arbitrary direction in phase space. The expression 



Q.^^^^ (6) 
^ trp2 trp2 ^ ^ 

is an intrinsic definition of phase space correlation, quite independent of the representation used to compute it and 
should not be confused with, e.g., local wave function correlations which do depend on the specific coordinate 
representation. 

In quantum optics it is customary to switch to the basis of creation and anihilation operators {q ± ip) / \/2h. In this 
context, the translation operator jSJ depends on the complex chords (^j, ± i£,q)/^/2h and is called the displacement 
operator 7, 8J. The semiclassical limit for a complex phase space is not as transparent as the real theory treated here. 
However it is quite feasible to effect phase space translations in the optical context 0- 

Some of the general formulae in the following sections are related to the theory of gaussian noise channels. Indeed, 
the "fidelity" of a gaussian channel is just the correlation © averaged with a gaussian window for the chords f |10| . 

III. THE CHORD FUNCTION AND THE WIGNER FUNCTION 

The chord symbol for an operator A is defined as 

A{0 = trf_^i , (7) 
allowing for the complete representation of A in terms of the unitary translations in phase space: 
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In the case of the density operator p, it is convenient to alter the normahzation, so that 

xiO - (^trf_,p (9) 

is the definition of the chord function, also known as the characteristic function, or the generating function. (In 
classical mechanics the displacement ^ results from a trajectory of which it is the chord - see Appendix fXl ) 
The Fourier transform of the chord symbol is the Weyl symbol A{X), 

^(^) = / TTTTTT e-'^''^^^A{0 , (10) 



in terms of the skew product 



XAC = P-^,-Q-^p. (11) 



In the case of the density operator, the Fourier transform of is the familiar Wigner function. But, since the Fourier 
transform of the translation operator T> itself corresponds to the classical reflection through the phase space point X 

El, 



.,XAi/n ^ 2^i?jf , (12) 



{2T:h) 



it follows that [TJl 



W{X) ^ j-l^trRxp . (13) 



Though translations and reflections are quite distinct operators, they combine to form the affine group in geometry 
|l3j |. which is transported into quantum mechanics by the operators Rx and |l2j| . The family resemblance is 
striking when viewed, for instance, from the position representation: 



2^J?x = / di 

whereas 

' dQ 



e»«p-QA . (15) 

It is well known that the Wigner function cannot be indentified with a probability distribution in phase space, even 
though this interpretation holds for marginal distributions and the calculation of averages of observables as phase 
space integrals. The main problem is that W{X) can assume negative values and indeed they are present for all pure 
states that are not gaussian (coherent) states. For this reason, it is also improper to refer to xiO as a charactersitic 
function, even though its derivatives do generate moments of q, p, and any polynomial function in phase space, taking 
proper care of the operator ordering; for instance, 

(9") = trq-p = {^nr ^ xiO ■ (16) 

For the case of the identity operator, we obtain the normalization condition: 

1 = trp = J dX W{X) = {2TThfx{0) ■ (17) 

In the case of operators representing observables with smooth, classical-like Wigner functions, e.g., polynomials in q 
and p, their corresponding Fourier transforms, i.e., the chord symbols are sharply localized (improper functions) close 
to the origin. On the other hand, the chord function for a normalized state is a proper function, and extends away 
from the origin, so as to represent truly quantum correlations. The process of decoherence, destroying the purity 
of the initial state, generically washes away the exterior structure of the chord function and compacts the mixed 
state onto the classical origin of chords. This has been shown for linear Markovian systems and will be the 
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subject of further work. Here we will describe the large and small scale features of the chord function and their close 
intertwining. 

Combining the definition H12|l of the Wigner function with the group properties of the translation and reflection 
operators jl2| it is easy to see that the Wigner function corresponding to the translated state ~ T^pf^^ is 

W^iX) = W{X-Tj) , (18) 

whereas the corresponding chord function is just 

XviO = e""^«/"x(0 . (19) 

Unlike the Wigner function, the chord function is not necessarily real, but it may be real for a particular choice of 
phase space origin. It is shown in Appendix IbI that the necessary and sufficient condition for this is that there exists 
a symmetry centre X, such that [p, Rx] — 0, and it is chosen as the origin. Since Rx has eigenvalues ±1, the pure or 
mixed state p, must then lie in the Hilbert subspace of either even or odd parity. For these parity-symmetric states, 
p±, the Wigner function and the chord function are obtained from each other by a mere rescaling: 

W^±(X) =± 2^X±(-2X) . (20) 

For general unsymmetric states, the real part of the chord functions is still determined by the diagonal part of p with 
respect to parity, whereas the imaginary part depends on the off-diagonal part. However, it is the intensity of the 
chord function, |x(^)P, that turns out to be most useful. 

The intrinsic definition of phase space correlations Cj in © is readily translated into the properties of Wigner 
functions and chord functions: 

tr pf^p fl = {2Trhf J dX WiX) W{X - C) ^ {2^n)^ J dr] e''^'^^'^ |x(r/)|' ■ (21) 

Thus the correlations of the Wigner function can be identified with and |x(77)P is just the power spectrum of 
W(x). 



IV. PURE STATES 



It is worthwhile to recollect some examples for which the chord function can be identified with the Wigner function, 
once the origin is translated to the symmetry center. In all the following cases we consider states of a harmonic 
oscillator with one degree of freedom and unit mass. 

i) Coherent states, jyy), are displacements of the ground state of the harmonic oscillator by Tr,. The Wigner function 
is just a gaussian centered on 77, 



Wr,{X) = ^exp 



(22) 



whereas, using lfTH|) and ifT^ . 



ir] 



exp 



1 

tlUJ 



(23) 



So, any translation of the coherent state merely alters the phase of the gausssian chord function that sits on the origin. 

(ii) A superposition of a pair of coherent states, j??) ± | — 77) is sometimes known as a "Schrodinger cat state". Its 
Wigner function is (here and below we set to = V) 



W±{X) 



1 



27r?i(l±e-'''/?i) 



h 



(24) 



It consists of two "classical" gaussians centred on ±77 and an interference pattern with a gaussian envelope centred 
on their midpoint. The frequency of this oscillation increases with the separation |2ry|. For the chord function. 



X±(0 



1 



47r?i(l±e-''V?i) 



+ e-(C/2+r,)V?. ± 2e-«'/4'' cos -^Ar, 

h 



(25) 
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this same configuration has to be reinterpreted. Now the internal correlations of the individual coherent states are both 
superimposed onto the neighbourhood of the origin, as in (i), while their cross-correlation generates new gaussians 
centred on the separation vectors ±2ry. Of course, the general case of coherent states |r?i) and |r?2) merely leads to 
gaussians centred on ±(771 — 772) with addition of the phase factor exp[i(77i + 772) A S,/2Ti]. 

(iii) Fock states, |n), i.e., the excited states of the harmonic oscillator, also have reflection symmetry with respect 
to the origin. Thus, from the exact Wigner function, first derived by Gronewold p^ . 

Wr.{X)^t}r,-x^/H(^^\ , (26) 



-fc t 

TT/i \ 

where L„ is a Laguerre polynomial, and H20() we obtain the chord function 

- (I) • 

It is interesting to note that the symmetry centre, which produces the maximum amplitude of the Wigner function 
is nowhere near the classical manifold with energy — + 5) However, this point lies in a region of narrow 
oscillations, so that it does not affect the average of smooth observables. 

All the above examples are singled out by some point of reflection symmetry, which must always be chosen as the 
origin for the chord function to be real. The chord function always assumes its maximum value l/(27r?i)^ at the 
origin, whatever the symmetry. For a pure state the proof is immediate because 

im^) = (2^;i)^x(-0 , (28) 

which cannot have modulus greater than one, whatever the symmetry. But even an average of overlaps cannot exceed 
one, so x(0) is also the maximum for mixed states. The Wigner intensity [VF(X)]^, need not have such a prominent 
peak in general. However we shall see in the section IVII that the large scale features of the semiclassical forms of 
the Wigner function and the chord function maintain a mutual correspondence, even in the absence of a reflection 
symmetry. 

General invariance with respect to Fourier transformation does hold for the correlation in the case of pure states. 
Indeed, combining (|28|l with H21(l we obtain 

{2^hrC, = |x(C)P = / e^'^-^/^Wiv)? ■ (29) 

This is a remarkable property of all pure states and is in no way restricted by special symmetry properties that relate 
certain Wigner functions to their respective chord functions. All the same, we can start by considering the example 
of the chord functions so far studied: (i) in the case of a single gaussian (I22|l the invariance is obvious, because the 
square modulus is a gaussian with the appropriate width for its Fourier transform to be of the same form, (ii) For a 
pair of coherent states, the chord function is a sum of gaussians. Its square modulus is also gaussian and we again 
return to the same function by Fourier tranformation. (iii) For Fock states this is not so obvious, but we can also 
write |x(OP ^.s 



(30) 



with coefficients that do not depend on So, the invariance of the ground state gaussian entails that of all Fock 
states. 

The Fourier invariance condition (|29|l includes as a special case the more familiar tracing over the full pure state 
condition — p. It can be easily checked that setting = in (|29|l gives the chord representation of the identity 
trp^ = tcp. It follows that the difference of both sides of H29I) for each chord f is a measure of the purity that 
generalizes the linear entropy. However, the loss of the phase information contained in the chord function, but absent 
in |x(^)|, implies that these are necessary conditions, whereas the full sufficient condition of purity is only — p, 
which is expressed in the chord representation as 

drixiv) X(C - Tl) e^^^"/"' - / drix^niv) x(C - v) = xiO , (31) 



with Xi/2{v) defined by l(T^ . 
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V. SMALL CHORDS 



We now consider the chord function and the correlation function of semiclassical states, i.e., those related to 
generalized Bohr-Sommerfeld energy levels. This section focuses on the limit of small chords. 
Starting from the Wigner-Weyl representation we may rewrite 



xiO - TT^z^ {f-d = I dX T^,{X)W{X) , (32) 



1 /f 
(27r?i)^ ^ {2TThy 

Besides the Wigner function, we have here introduced the Weyl symbol for the translation operator: 

2^ tr f _^ Rx = T„^ {X) = e-'^''^/'' . (33) 

If ^ is small enough, i.e., |^| < h, then T^(X) behaves like a smooth, classical-like symbol. It still oscillates, but with 
a classical wavelength. In such a case, to a good approximation, we can replace the Wigner function in H32|l by the 
simplest semiclassical expression 

Wj{X)c,-^S{I{X)~I) . (34) 

Here I(X) is the set of L action variables for an integrable system with L degrees of freedom and T is the set 
of quantized action values for this particular state [l^ Il9j | . Within this approximation, the average of a quantum 
observable A is just a purely classical average over a torus: 

where 9 are the angle variables conjugate to the actions, which describe positions on the quantized torus, and A{9) = 
A{X{6)). In the case of the representation for the chord function H32f) . we obtain 

- {2nh)L J (27r)^ ■ ^ ' 

Certainly, a bad choice of origin will lead to large phases in (|36|l , but we have already studied the trivial phase change 
due to translating the origin. Thus, we can increase the quality of H3()l) by choosing the origin to minimize |^(0)| on 
average. 

Let us check the semiclassical approximation (|36ll for the simplest case of the Fock states, discussed in the previous 
section. Choosing w = 1, we have action- angle variables that are merely canonical polar coordinates 

(7 = \/2X cos 6* , Vol sine. (37) 

Thus, choosing ^ along the p-axis, without loss of generality, we obtain 

^(^^ " 2^/ S e-^l«l-^«/^ = —Jo I 1 , (38) 

where Jq is a Bessel function. This result should be a good approximation to the exact formula l|27|l for |^| ^ fi and 
large n = 2/fi — 1/2. To check this, first note that x{C) oscillates with a wavelength A ^ h/^/T. In terms of this scale, 
and defining z = I^JI/A, the argument of the Laguerre polynomial in (|27|l reads 

Hence, for <C we recover (|38f) by using the formula ^3 

In the next section we shall study the semiclassical limit of the chord function for large chords and show that in 
the particular case of the harmonic oscillator this coincides with the expansion of H38|) for large argument: 

Jo(.)--|=cos(,-J) . (41) 
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VI. SEMICLASSICAL THEORY 



The easiest path to obtain the semiclassical form of the chord function for finite chords is to start from the general 
WKB wave functions. The underlying classical structure is assumed to be a curve (for L — 1) or a Lagrangian surface 
{L > 1) defined by the actions /. Fixing the action, such structure can be described locally as a function p = p{q, I)) 
which may have several branches. Thus we define 



p{q',I)dq' 



which is the generating function for the transformation {p, q) — > (/, 9), such that 

dl ' dq ^ 
Then the corresponding WKB wave functions are linear combinations of 



iili'i 



dot 



d'S{q,I) 



dqdl 



1/2 



^lS{q.I)/h 



(42) 



(43) 



(44) 



for the various branches of S{q,I), with c a normalization constant [TM Il9||. 

The similarity between equations (|14|) and H15|l for the operators that determine the chord function and the Wigner 
function, allows us to follow the same steps as Berry (16ii] for the semiclassical Wigner function. Thus we presume 
that 



s[Q' + ^,i)-s(q'-^,i]-^,q' 



1/2 



(45) 



dqdl dqdl 
has a stationary point within a single of branch of S{q, I). The stationary phase condition for Q' — Q 

Px(q + |)-Pi(q-|) (46) 

specifies that the geometrical chord corresponding to the arc of pi{q) = p{q, I — I), lying between q = Q — S.q/'^ and 
q = Q + coincides with the given chord ^. 

Unlike the Wigner function, the phase of the chord function depends on the choice of phase space origin, but this 
phase is trivially specified by the overall phase factor (|19|l for a translation by a vector rj. Thus, let us for now assume 
that the origin is translated to (0,(5). Then the stationary phase action is just the area shown in Fig. ^a). The 
alternative constructions (b) and (c), also shown in Fig. ^ have the same area and enjoy the advantage that they 
survive arbitrary linear canonical transformations which preserve the origin. In other words, we may always consider 
the phase for the semiclassical chord functions to be defined by a sum of areas: (i) the area sandwiched between the 
chord and the arc into which if fits, S{X,T). This is just the same as the one for the semiclassical Wigner function, 
except that there the construction starts from the centre of the chord X [Fig. H^b)] , rather than from the chord itself, 
(ii) The area that is added to this may be taken as that of the parallelogram obtained by tranporting the centre of ^ 
to the origin [Fig. Q^b)] , or the polygonal figure of eight constructed by the chord and its reflection around the origin 
[Fig.njc)]. In either case the value of the area is just X A ^. Therefore we construct the chord generating function 



S{tI)^S{X,I)~XA^ 



(47) 



which is just the Legendre transformation of the centre action, appropriate to the Wigner function, as discussed in 
Appendix El This geometry is immediately generalized for L > 1: the action for any polygonal figure is just the 
algebraic sum of the areas of its projections on each conjugate plane. Furthermore, it does not matter which arc is 
chosen between the tips of ^ along pxiq) because this is a Lagrangian surface. 

Thus, the contribution of each realization of the chord ^, with centres Xj and closing arcs 7j on the quantized curve 
or Lagrange manifold I{x) = X has the form 



Xj(e) =^j(Oexp 



(48) 




FIG. 1: The stationary phase condition for the semiclassical chord function xiO is just that ^ be a geometrical chord for the 
curve I{x) = I, locally given by px{q). The stationary phase itself is given by the shaded area in (a). The constuctions in (b) 
and (c) are invariant with respect to linear transformations and have the same areas. 



Here aj is the signature of the matrix 



[S{q+,I)-S{q-,I)] , 



(49) 



where q± = Q ± ^/2, so that in the case of the curve in Fig. ^we have cr = +1. 

The amplitude for each reahzation of the chord, within a normalization factor, is given by 



det 
det 



dpi 
dq 



dpi 



dl dl 
det — {p{q+),q+)) det — (p(g-), q-)) 



dl 


dl 


dl 


dl 






dp 


^ dq 


dp 


^ dq 







because 



dl_ 

dq 



dq 



dl dp 
dp dq 



= 



(50) 



(51) 



along the classical surface. These amplitudes coincide with those of the Wigner function evaluated at the point 
Xj{^). In the case where L = 1, the 1x1 determinant is only a single factor. Taking I{x) to be a Hamiltonian, the 
corresponding phase space velocity tangent to the phase space curve is 



±1 — J 



m 

dx 



where J is the symplectic matrix (|A4ll , and we may interpret the amplitude of the chord function as 

\MO\-'^{ii)+^{xi)_ = {i+,i-} , 

where the displaced actions are defined by 



nx±l 



(52) 



(53) 



(54) 



(see Fig. El. In the general case each element of the determinant is such a Poisson bracket for the different L action 
variables |20| . 

For a maximal chord, termed a diameter, such that i;_ becomes parallel to x^, the amplitude diverges to generate 
a caustic of the chord function. Beyond this boundary x(0 becomes negligible. In the limit of small chords, x± 
are almost parallel, and the amplitude may be arbitrarily large. In this region the above semiclassical theory ceases 



FIG. 2: The semiclassical chord function, xiO^ is determined by the geometrical chord ^, centered on X, and joining the points 
X- and x+ on the quantized classical curve. The amplitude depends on both phase space velocities at the tips, as explained in 
the text. 



to operate and, because the chord normalization condition is precisely x(0) = l/27r?i, we have neglected to define 
the overall normalization factor. Actually this limit around the origin is much more nasty than normal semiclassical 
caustics, since the entire classical manifold can be defined by a succession of infinitesimal chords. It is thus a nontrivial 
problem to construct a semiclassical uniform theory including large and small chords. 

However, it is possible to connect the semiclassical theory to that for the small chords of the previous section. Let 
us take the case that the given chord, ^, is parallel to the q axis, so that we can rewrite (|45|) as 



X(0 



using H43|l . Expanding 



s[q + ^,i)-s(q-^,i 



, 961 



1/2 



, de 



1/2 



(55) 



(56) 



we obtain 



, de 



1/2 



, de 



1/2 



, de 



[i+o[e)] 



(57) 



whereas 



S[Q 



(58) 



Hence, we can cancel the amplitudes in (I55|l by changing the integration variable, Q — > 0: 

1 



X(0 



{2TTh) 



2-K in 



L I, 27T 



(59) 
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The general form of (|59|l holds for any choice of origin and direction of the small chord ^. 

We have here rederived the small chord expression for the chord function used in the last section in a way that has 
several advantages. First we notice that Berry's original derivation of the 5-function approximation for the Wigner 
function also linearized the action for small chords 16]. But here the chord is indeed fixed and can be assumed to be 
small, whereas there the integral is over all chords. We should note that (|34() . or the extension of H59|l for all chords, 
is tantamount to approximating 

P^2-l (60) 

We can now evaluate the error in (|59|l when the chord ^ is large enough for both this integral and H45|) to be evaluated 
by stationary phase. Consider the convex curve I{x) = X in Fig. |3fa). The difference in the phases for the stationary 
phase evaluation at each stationary point of (|59|l or the full integral 145|l is the part of the area near the corner of 
the circunscribed parallelogram with two sides ^ tangent to the curve, lying outside the curve. If we approximate the 
curve as a parabola around the point of tangency 5p = a{Sq)^, the leftover area is just ^(<^) — 2a|^p/3. So, if the 
diameter between the parallel tangents with the direction ^ is C{0> true stationary area for each realization of the 
chord is just 

^(0 = ^CA^-^. (61) 

Thus the stationary phase evaluation of the full integral is valid as long as |C A ^| > h, whereas the integral (|59|l can 
be used if a|^p < Ti. Therefore the range of overlap for both approximations is 

^<m<(-Y\ (62) 



ICI 



so that we can always make the transition between both approximation as ?i — > 0. 



VII. CONJUGATE CHORDS 



The pair of parallel tangents shown in Fig. 3(a) can be defined as the realizations of an infinitesimal chord for the 
curve I{x) = X. The points where the tangents touch this curve are joined by the diameter In the case of a finite 
chord ^, we can also define the conjugate chord, r/, as that which closes off the inscribed parallelogram with the sides 
^, as shown in Fig. Olb). Clearly, the area of this parallelogram is ^ A 77, which is smaller than ^ A C, used in the 
previous section. Choosing the origin at the centre of this parallelogram, we obtain the phase of each contribution to 
xiO as f A 77/2, added to each of the chord areas for ^. 

Increasing the length of ^ while keeping its direction constant, the pair of realizations eventually coalesce along a 
diameter of the closed curve. Thus we obtain in this limit an interchange in the role of the "conjugate chords" and 
77. Now it is the small chord 77 that determines the phase of the contribution as the chord ^ approaches the caustic of 
xiO its diameter. 

It is not only in the geometry of the curve that the conjugate chords ^ and 7/ interact. Indeed, the phase difference 
between both realizations of the chord ^ coincides with the sum of the chord areas for the realizations of the chord 
77. To see this, notice that this phase difference is just the shaded area in Fig. EJa) divided by Planck's constant. 
But, because this curve is quantized (its enclosed area is an integer factor of 27r?i, plus a Maslov correction), the 
phase difference for is the same as for the unshaded area. Notice now that the phase difference can also be 
obtained by translating the whole quantized curve by the chord ^ and measuring the overlap. This is just the classical 
correspondence for the definition of x(C) of a pure state in H28|) . If the classical quantized curves all have a symmetry 
centre, it makes no difference whether or not the curve is reflected about its centre prior to translation. In this way we 
obtain a semiclassical version of the trivial relation between the Wigner function and the chord function in the case 
of centre-symmetrical systems. The overlap between the quantized curves, or quantized tori, was already employed 
to study the semiclassical approximation to the Wigner function [20|. 

Both the Wigner functions and the chord function can be considered as special cases of the overlap of two semiclas- 
sical states lip) and |(/)), corresponding to classical curves, or Lagrangian surfaces. In the case of the Wigner function 
1^) = Rx\ip), whereas, for the chord function |(/)) = T^\i{j). In the general formalism of Littlejohn each inter- 
section of the manifolds corresponding to and |0) determines a semiclassical contribution. The phase difference 
between these contributions is determinated by the overlap area, just as in Fig.^Jb), within a Maslov correction. The 
amplitude of each semiclassical term is then given by Poisson brackets between both manifolds at the intersections. 
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FIG. 3: (a) The difference in the pair of stationary phases for Ij59|l is proportional to the paraUelogram that has ^ as one of its 
sides and is doubly tangent to the quantized curve. The difference between this area and the pair of stationary areas in the full 
semiclassical formula is just the shaded area at the corners, (b) A pair of stationary phase realizations in a convex quantized 
( 
( 




FIG. 4: The phase difference for the pair of realizations of the chord ^ is the shaded area in (a) divided by h. The quantization 
of the curve equates this phase (modulo 2n) to the phase for the shaded area in (b), obtained by the overlap of the curve 
with its translation by ^. Adding to this the area of the inscribed parallelogram, we obtain the phase difference for the pair of 
realizations of the chord rj in (c). 
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Obviously, in the present case, this coincides with the Poisson brackets for the same manifold at either tip of the 
chord, as obtained in (|53|l . Perharps this general point of view provides the more pleasing explanation for the identical 
form of the Wigner and the chord amplitudes. 

The relation between conjugate chords, ^ and 77, and the corresponding evaluations of the semiclassical chord 
function for the same state j?/'), i.e., XV'(0 ^^'^ Xtpiv)j is mediated by the parallelogram formed by ^ and rj. The 
phase difference AS'(^), in units of h, between both contributions to XipiO can be pictured as either the shaded area 
in Fig.^a), or the unshaded area corresponding to both realization of the chord rj. But, if we now add ^ A 77 to these 
unshaded areas, we obtain the new shaded area AS{ri) in Fig. ^c), which determines the phase difference of both 
contributions to Xi'iv)- 

Squaring the semiclassical approximation for xiv) '^c obtain 

{27ThfC, = \xivf ^ {I+,I-K' +{I+,I-}^' + {I+,I-}r{I+,I-}2' cos (63) 

in the simple case where I{x) — 2 is a convex curve, so that there is a pair of realizations for the chord, 77, with 
intensities given by (|53() . Inserting H63|l this into the Fourier integral (|29|l . we find that the first two smooth terms 
only contribute to the classical neighbourhood of the origin. It turns out that the stationary condition for 77 in (|29|l for 
a given value of ^ in the last term is precisely that 77 be the conjugate chord to ^, i.e., that ^ and rj form an inscribed 
parallelogram in the curve corresponding to {i/j). This is a consequence of equation (|46() . Therefore the Fourier 
invariance of the quantum correlation results semiclassically from the relation of the chord function xiO ^ xMi for 
conjugate chords. 

In the limit as the family of chords with a fixed direction approaches its maximum value for a given quantized 
curve, i.e., a diameter, the conjugate chord approaches the origin in the direction of the parallel tangents at the tips 
of the diameter. However, in this region, the simple semiclassical amplitudes Aj given by H53|) become singular and 
each of the three terms in (|t)3|l contribute. It is here necessary to replace the semiclassical contribution by uniform 
approximations in terms of Airy functions. This will be the subject of further work. 

The role of diameters as conjugate chords near the origin is brought forth by combining H29|) with H36() . for |^| — > 0: 

^ J {2tt)L {2tt)L ^ ' 

We can interpret x{9+) — x{9-) as the set of chords supported by the quantized curve, or L-dimensional torus. This 
confirms the role of the first two terms of H()3|) which contribute to the Fourier transform of for small ^: all the 
chords supported by the curve contribute. However, as soon as ^ grows in modulus enough that we may evaluate 
by stationary phase [condition H62(l ]. the dominant contribution comes from the diameter C(0; i-^-; C — 
such that 

in other words, the tangents at 6± are parallel. 

Up to now we have analyzed cases where the classical structure is a continuous curve or surface in phase space. Let 
us now consider an alternative classical setting, a state 

j 

where |?7j) are coherent states centered on the phase space points rjj. Then the density operator is 

and its chord representation Xi'iC) is a simple generalization of H24|l : each chord ±(r/j — rjk) is the centre of a gaussian, 
whereas the diagonal terms in 1671) interfere collectively in the neighbourhood of the origin. Taking |Xj/'(OI^: again 
have gaussians centred at rjj — rjk and at the origin. It is easy to see that the Fourier transform merely interchanges 
the contributions to the origin with the pair of gaussians centred at ±(?7j — rjk). However, there is a new contribution 
to the Fourier transform if the four vectors 771, 772, 773, 7?4 form a parallelogram, i.e., 772 — 771 ~ 773 — 774 and consequently 
773 — 7/2 — 774 — 771, where ~ means that the vectors differ by 0{\/fi). Then the Fourier transform of |x(C)P for 
^ = 771 — 7/4 receives contributions from the neighbourhood of 772 — 771 and vice- versa. Thus, again we verify the role of 
conjugate chords in the Fourier invariance of the quantum correlations. The consistency of fitting semiclassical states 
with gaussian coherent states has been recently investigated in 22] . 
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VIII. DISCUSSION: RESURGENCE OF PURE STATE CORRELATIONS 



We recapitulate our findings: As we displace a pure state that is classically extended, i.e. that quantizes with a 
large quantum number, the correlation in a given direction goes through four different stages: a) An initial, relatively 
simple stage of very short chords in which the behaviour is quadratic in the displacement and determined purely by the 
phase space extent of the state. This is the region that has received most attention in the litter ature so far pll5ll23l|. 
b)A second oscillatory stage ruled by the points of intersection of the two displaced tori, with relative phases that are 
semiclassically determined, c) A third stage characterized by a chord caustic where the semiclassical contributions 
diverge and where uniform approximations are still needed, d) An asymptotic region where the correlation decays 
uniformly to zero. We have here given special attention to the transitions between stages (a) and (b) . Other large 
systems such as a widespread superposition of coherent states will also have large classical values for (p^) or (q^). For 
all such pure states the initial decay of the correlations with growing displacement will be followed by their oscillatory 
resurgence for those chords that cause the underlying classical structures to overlap. 

It is important to consider higher phase space dimensions, since all our examples were restricted to a single degree 
of freedom. Again, it is clear that in the case of arbitrary superpositions of pairs of coherent states the straight Fourier 
analysis of the Wigner function leads to very similar pictures for the pair of classical-like gaussians in the Wigner 
function and the peaks of correlations at the chords that separate them. 

Integrable systems of higher dimension might appear to be harder to analyze, but this is not so. It was show in 
Ref. |2flj that the Wigner function corresponding to L-dimensional quantized tori are characterized by a Wigner caustic 
of dimension 2L — 1 in which the torus appears as a higher singularity. At this boundary between an oscillatory inner 
region of phase space and the evanescent region outside, the Wigner function attains an amplitude maximum, which 
can be described locally by an Airy funtion and its derivative. In the case that there exists a reflection symmetry 
centre, the chord function must be identical to the Wigner function within a phase and the rescaling (|20|l . Even 
without such a reflection symmetry, the region where a finite L-torus and its rigid translation intersect is bounded by 
a {2L — l)-dimensional surface at which they touch nontransversally - the chord caustic. 

Finally we must consider the case of an eigenstate of a chaotic Hamiltonian. In this case Shnirelman's theorem 
guarantees that most states are ergodic in the sense that the average of smooth functions of the observables p and q 
are given by a classical average over the energy shell. In its simplest form, quantum ergodicity may be taken as the 
Berry- Voros hypothesis that the Wigner function is approximately a Dirac delta function on the energy shell '255 ■ 
It follows immediately that for small chords we may adapt the discussion in section V to obtain the chord function 
for the n'th energy eigenstate as 

Xn{0-{e )n-J^ Jdx6iH{x)-E,,) 

(see also 0,113). Conjugate to this short chord behaviour there must be appreciable large arguments for which Xn{£,) 
is also large, corresponding to the largest chord fitting between a pair of points in the energy shell for any given phase 
space direction. Once again, in the case of a refiection symmetric chaotic Hamiltonian, we may invoke the (scaled) 
identity with the Wigner function which is dominated by the energy shell itself, according to Shnirelman's theorem 
0. 

In all the extensions to higher dimensions, large chords ^, generate oscillations in the Wigner function with wave 
vector J^. In the case of a semiclassical state constructed on a Lagrangian torus, the fine structure in the Wigner 
function results from the interference of a finite number of locally plane waves. The fine structure of ergodic chaotic 
Wigner functions has not yet been explored. 

It should be noted that our general description is in no way limited to stationary states. The curves and surfaces 
which we have mainly treated may be evolving classically, while the corresponding quantum system also evolves. We 
have not treated here this dynamics, but our description is valid for any snapshot. Even for the evolution of a localized 
wave packet in a chaotic system, it is possible to advance that long scale quantum correlations develop as the classical 
packet spreads over the energy shell. 

This study of the coherence properties of pure states suggests that it is not always profitable to transform back to the 
centre phase space of the Wigner function from the phase space of chords, since it is here that the full structure of chord 
conjugacies is manifest. Furthermore, there is a decided advantage to allow the collapse of all classical information 
onto the neighborhood of ^ = 0. All that extends out in the chord phase space are signs of quantum coherence and 
these are structures that are lost in the nonunitary evolution that results from tracing out the interaction with an 
uncontrollable environment. 
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APPENDIX A: CONJUGATE PHASE SPACES 



The points in phase space for a system with L degrees of freedom are here denoted simply by a; = (j), q) = 
(pi, • • • ,pl, gi, • • • , (7l). In the classical limit, ?i ^ 0, we may associate a quantum state to such a phase point, but we 
need pairs of points (a:_, to describe operators A, such as or {p+\A\p-). In both these cases, we only 

use, in fact, half of the phase space variables, as consistent with the uncertainty principle. 

Another alternative is to use either the centre 

X = {P,Q)^{x++x-)/2 (Al) 

or the chord 



^={^p,^g)^X+-X^ . (A2) 

As shown in 0, we may then identify the argument in the Weyl symbol of A, A{X), with the center IjAip . 
In the same way, the argument in the chord symbol A(^) corresponds to the chord l|A2|l . In the case of a unitary 
transformation, corresponding classically to a trajectory, we identify ^ as the chord corresponding to the arc {x{t), < 
T <t} joining x_ = a;(0) and x+ = x{t). Of course, X is the reflection centre for this pair of phase space points (12| . 

Canonical transformations in phase space are generated implicitly by the centre action, i.e., the generating function 
S{X) such that 

.Ji = §, (A3) 

where 



or by the chord action S{^), such that 



-JX=|. (A5) 

Comparing with more familiar generating functions, e.g., S{q-,q^), 

dS dS 

it follows that the conjugate variable to X should really be rather than the chord itself. However, this lacks a clear 
geometrical interpretation in terms of x±, so it is better to keep ^, but replace all scalar products by skew products. 
Hence, in all Fourier transforms we use 

L 



APPENDIX B: PARITY EIGENSTATES 



Any state \ip) in Hilbert space can be decomposed into components of even or odd parity (eigenvalue +1, or —1) of 
any of the reflection operators Rx- Indeed these projectors were presented by Royer ^J] as 

= i (l ± TThRx) , (Bl) 
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so that the parity decomposition of an arbitrary density operator, p, is 

pX c. pX 

tr p 

Since these reduced density operators commute with the reflection operator, the Wigner function corresponding to 
RxP is 

W^{X') = 2HvRx'RxpI = ± 2^tri?x'pf = 2^6^' '''^''/''trfafx'-x)^^ , (B3) 

where we used the general group relations between translations and reflections ,12] , Thus, shifting the origin to the 
centre of symmetry, we obtain 

W^{X) = ±2^X^{-m- (B4) 



In the case that is the projection of an arbitrary p according to (jB2p . then the general form of the projected 
Wigner function given by [26| is transported to the chord function as 

± x(0+x(-0±2W^(g/2) 

^° ^ A[i±.nwm — ■ 

Evidently, all chord functions with pure parity are real. The reciprocal is also true, because the imaginary part of 
x{Ct cancels if and only if 

tr(f_^p) = x{0 = X(e)* = = trf^/S • (B6) 

But, using the group properties of translations and reflections, 

trf_jp = trRof^Rop = trf^RopRo , (B7) 

which is only equal to (|B6I) for all ^ if [p, Rq] — 0. 

A mixture of pure states, each of which has definite parity, commutes with Rq and hence produces a real chord 
function. The imaginary part of the chord function is related to the off-diagonal parity representation of the density 
matrix, i.e., for an orthogonal basis of odd and even states, this is the block of the density matrix coupling the different 
parities. 
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